Polarization of prompt J/V' in PP J/ip + X at ^/s = 200 GeV 
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Within the framework of the nonrelativistic QCD factorization approach, we compute the cross 
section and polarization of prompt J/ip produced from proton-proton collisions at the center-of- 
momentum energy ^/s = 200 GeV. We present the transverse-momentum distribution in the forward- 
rapidity region 1.2 < \y\ < 2.2 and the rapidity distribution over the transverse-momentum range 
2 GeV < pT < 20 GeV. The perturbative contributions are computed at leading order in the strong 
couphng constant. We predict sUght transverse polarization of J/ip in the forward-rapidity region, 
while that for the midrapidity region is shghtly longitudinal. The transverse-momentum distribution 
agrees well with the PHENIX preliminary data and the color-singlet-model prediction at next-to- 
leading order in as, but disagrees with the result from the leading-order color-singlet model or the 
s-channel-cut method. 

PACS numbers: 12.38.-t, 13.85.Ni, 13.88.-|-e, 14.40.Pq 



I. INTRODUCTION 

As one of the long-standing puzzles in heavy- 
quarkonium phenomenology, the polarization of prompt 
J/ip produced with large transverse momentum (pt) at 
hadron colliders has not been clearly understood yet. Ac- 
cording to the nonrelativistic quantum chromodynamics 
(NRQCD) factorization approach [1], the dominant pro- 
duction mechanism of J ftp with large px in hadron colli- 
sions is the gluon fragmentation into the color-octet spin- 
triplet 5- wave heavy-quark-antiquark pair [QQsi'^ Si)] 
that evolves into the meson [2] . Because the correspond- 
ing long-distance transition preserves the spin at leading 
order in v, which is the typical velocity of the heavy 
quark of the meson, the J/'ijj should be transverse at 
large pt [3]. However, the run II measurement by the 
Collider Detector at Fermilab (CDF) Collaboration [4] 
agrees with neither the run I data [5] nor the NRQCD 
predictions [6-10], although the feeddowns from the P- 
wave spin-triplet states Xcj for J — 0, 1, and 2 dilute 
the transverse polarization [9, 10]. The problem has not 
been resolved even after considering corrections in vari- 
ous aspects [11-22]. Therefore, it is worthwhile to see if 
NRQCD still confronts with experiments of lower center- 
of-momentum (CM) energies ^/s, where the dominance 
of the gluon fragmentation may not happen. 

In this aspect, the inclusive J/ip production in pp colli- 
sions at y/s = 200 GeV, where the fragmentation does not 
dominate [23] , at the Brookhaven Relativistic Heavy- Ion 
Collider (RHIC) may provide independent experimental 
constraints to the J/ip polarization. Because only a tiny 
fraction of the PHENIX inclusive J/ijj samples contain 
the J/tp from the B decay [24], we can treat this as 
prompt J/ip- In the case of the cross section, the NRQCD 
predictions in Refs. [23, 25] are consistent with the run-3 
PHENIX data [26, 27]. As a new production mechanism, 
the s-channel-cut approach [28] , which suffers from a crit- 
icism [29], also explains the data [30]. The color-singlet- 



model (CSM) prediction for the direct J/ip production 
rate has been computed at next-to-leading order (NLO) 
in the strong coupling [31, 32], including contributions 
from eg fusion. The NLO CSM prediction agrees with the 
data for the rapidity distribution, assuming that about 
40 % of the J/ip samples are from higher resonances, 
while it underestimates the data for the pt distribution. 
The NLO NRQCD prediction has been found to agree 
with the data for the pT distribution [20, 21]. In addi- 
tion to the cross section [33], the PHENIX Collaboration 
has measured the polarization [34] of the inclusive J/tp. 
In the midrapidity region \y\ < 0.35, the PHENIX data 
for the J/ip polarization [34] agree with the predictions 
of both NRQCD [35] and the s-channel-cut method [36]. 
However, the PHENIX preliminary data in Ref. [37] for 
the forward-rapidity region 1.2 < \y\ < 2.2 disfavor the 
s-channel-cut prediction [36, 37]. The PHENIX Collabo- 
ration is currently carrying out a comprehensive analysis 
with the data from the forward-rapidity region. There- 
fore, it is desirable to provide NRQCD predictions that 
can be compared with the forthcoming updated PHENIX 
data for the inclusive J/ip polarization in the forward- 
rapidity region. 

In this paper, we present an updated NRQCD pre- 
diction for the polarization of prompt J/tp produced in 
the forward-rapidity region 1.2< \y\ < 2.2 in pp colli- 
sions at y/s = 200 GeV. The perturbative contributions 
are computed at LO in as, like the previous analysis in 
Ref. [35], in which the cross section and polarization in 
the midrapidity region \y\ < 0.35 within the px range 
1.5 GeV < Pt < 5 GeV are predicted in detail. The 
PHENIX Collaboration is also preparing to release the 
rapidity dependence of the J/ip in both the midrapidity 
and forward-rapidity regions. In this ongoing analysis, 
they are using a conservative lower px cut pt > 2 GeV 
to avoid the region where LO NRQCD may break down. 
We also present a NRQCD prediction for the rapidity 
dependence of the polarization of prompt J/^ in the re- 
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FIG. 1: The difTerential cross section of prompt J/'i/) with the 
rapidity 1.2 < \y\ < 2.2 in pp colhsions at y/s = 200 GeV as 
a function of pr- The shaded band represents the NRQCD 
prediction and the band surrounded by a sohd curve is the 
color-singlet contribution at LO in Oa- 

gion \y\ < 3. Except for these changes in the kinematic 
region, our analysis presented here is very similar to that 
in Ref. [35] so we do not reproduce details of our strategy 
to compute observables and explanations on the input 
parameters given in Ref. [35] . 

The remainder of this paper is organized as follows. 
We first present the NRQCD prediction for the pr de- 
pendence of the differential cross section followed by the 
NRQCD prediction for the pr dependence of the polar- 
ization parameter a integrated over the forward-rapidity 
region in Sec. II. The rapidity distributions of the cross 
section and polarization parameter a integrated over the 
Pt range 2 GeV < pr < 20 GeV are given in Sec. Ill, and 
we summarize in Sec. IV. 



II. PT DISTRIBUTIONS 

In this section, we provide the NRQCD predictions for 
the Pt distributions of the production rate and polar- 
ization parameter a for the prompt J/'ip produced in pp 
collisions at -^/s = 200 GeV. The basic strategy to com- 
pute the cross section and polarization parameter a can 
be found in Ref. [35]. 

The cross section for the prompt J/tp in pp colli- 
sions is computed by including the LO parton processes 
ij — cc + k, with i, j — g, q, q and q ^ u, d, s. We 
use numerical values for the NRQCD matrix elements 
(ME) given in Ref. [9] . For the parton distribution func- 
tions (PDF), we choose MRST98LO [38] as the default 
value and CTEQ5L [39] for comparison. The transverse 
mass rriT — (4r7i^ -t-pf^)"'^/^ is used for both the factoriza- 
tion and renormalization scales ^ with the charm-quark 
mass rric = 1.5 GeV. is evaluated from the one- loop 




-0,2 r 
-0.4 r 
-0,6 - 

1,2<lyl<2.2 

-0.8 - 



12 3 4 5 

Pt (GeV) 

FIG. 2: The polarization parameter a for the prompt J/tp 
with the rapidity 1.2 < \y\ < 2.2 in pp colhsions at — 
200 GeV as a function of pr- 



formula using the value of Aqcd given in each PDF 
set [38, 39]. To estimate theoretical uncertainties, we 
follow the method given in Ref. [35] . 

Our prediction for the differential cross section of 
prompt J/ip is shown in Fig. 1 as a function oipT- The 
shaded band indicates the LO NRQCD prediction and 
the LO GSM contribution is displayed as a band sur- 
rounded by a solid curve. The rate is averaged over the 
forward-rapidity range 1.2 < jyj < 2.2, and its explicit 
normalization is given by 

-^^{d?a/dydpT)y = [ f' f dy, 

2TipT 2ttpt X 2 7i.2<|y|<2.2 dydpT 

where Bee is the branching fraction for J /ijj — > e+e^ and 
the factor 2 in the denominator on the right side is from 
the range of integration 1.2 < \y\ < 2.2. 

The result shows that more than 70 % of the prompt 
J/i/j's are direct: At px = 1.5 GeV (5 GeV), the contribu- 
tions from direct J/V', feeddown from ip(2S), and feed- 
down from XcJ are 83% (71%), 8% (11%), and 10% 
(19%) of the prompt J/tp^s, respectively. According to 
Fig. 1, the LO color-singlet contribution is negligible, 
which amounts to 4% (2%) of the direct J/ip, while the 
color-octet ^Pj or '^Sq and channels are 94% (81 %) 
and 1.5% (17%), respectively, at pr = 1.5 GeV (5 GeV). 
The uncertainties are evaluated by considering the vari- 
ations of nic = 1.50 ± 0.05GeV, rriT/^ < H < 2mT, 
NRQCD ME, and PDF [38, 39] as described in Ref. [35]. 
The uncertainties from /i of 83 % (98 %) dominate over 
those from nie [10 % (2 %)] and those from the NRQCD 
ME and PDF [7% (0%)] at pr = 1.5 GeV (5 GeV). 

The polarization parameter a is defined by a = (ctt — 
2(Tl)/{(Tt + 2(Tl), where (Tt (ctl) is the cross section for 
the transversely (longitudinally) polarized J/V'- If J li^ is 
completely transverse (longitudinal), then a = +1 (— 1). 
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FIG. 3: The diflferential cross section of prompt J/ij) with the ^IG. 4: The polarization parameter a for the prompt J/t/; 
transverse momentum 2 GeV < pt < 20 GeV in pp colhsions ^j^h the transverse momentum 2 GeV < pr < 20 GeV in pp 
at V5 = 200 GeV as a function of y. coUisions at = 200 GeV as a function of y. 



If J/V' unpolarized, then a — Q. The polarized cross 
sections depend on the frame. We choose the hadron 
CM frame that was employed in the PHENIX analysis. 

The pt dependence of the polarization parameter a in- 
tegrated over the forward-rapidity region 1.2 < \y\ < 2.2 
is shown in Fig. 2 with the same style as Fig. 1. The 
NRQCD prediction is slightly transverse over the whole 
range of px, while the color-singlet contribution gets 
more transverse as pt increases. In comparison with the 
prediction for the midrapidity region, the uncertainties 
mostly (> 98%) come from the NRQCD ME and PDF. 

Our results for the polarization agree with the 
PHENIX preliminary data [37] and are also compati- 
ble with the CSM prediction at NLO [32], which has 
larger uncertainties than the NRQCD prediction. The s- 
channel-cut prediction disagrees with both NRQCD and 
the PHENIX preliminary data [37]. 



III. RAPIDITY DISTRIBUTIONS 

In this section, we provide the NRQCD predictions for 
the rapidity distributions of the production rate and po- 
larization parameter a. The PHENIX analysis of the 
rapidity dependence of the J/ip is under way with a con- 
servative cut Pt > 2 GeV which avoids the kinematic re- 
gion where the fixed-order calculations may break down. 
Because the differential cross section at px = 5 (20) GeV 
is less than 10 % (0.1 %) of that for pr ^ 2 GeV, the to- 
tal cross section should be insensitive to an upper limit 
> 10 GeV and the neglect of the gluon fragmentation 
can be justified. Our results show that about 79 % of the 
total production rate comes from the direct J/ip, and 
feeddowns from Xcj and i/'(2S') occupy 13% and 8%, 
respectively. In Fig. 3, we show the NRQCD prediction 
for the differential cross section integrated over the range 



2 GeV < Pt < 20 GeV as a function of y. According to 
this figure, the LO color-singlet contribution is negligi- 
ble over the rapidity range \y\ < 3. The lower cut for 
Pt severely eliminates the color-octet and color-singlet 
■^5*1 contributions to occupy only a few percent of the to- 
tal production rate, making ^Pj and ^Sq octet channels 
dominate. The dominant source of the uncertainties in 
the rapidity distribution is the scale dependence which 
occupies about 92% (96%) at y = {\y\ = 3). The con- 
tribution of TOc dependence is about 6 % (3 %) and that 
of NRQCD ME and PDF is about 3% (1%) at y = 

(M=3). 

The rapidity distribution of a integrated over the range 
2 GeV < PT < 20 GeV is shown in Fig. 4. The LO CSM 
prediction is always transverse over the whole range of 
the rapidity. The central values of the NRQCD pre- 
diction are slightly longitudinal at midrapidity region, 
and there is a turnover around |?/| ~ 1 from which the 
NRQCD prediction becomes transverse. The uncertain- 
ties of the NRQCD prediction become maximum (mini- 
mum) at y = (|y| w 1). Therefore, the uncertainties of 
the NRQCD prediction for the forward-rapidity region in 
Fig. 2 are significantly smaller than those for the midra- 
pidity region presented in Fig. 2 of Ref. [35]. More than 
96% of the uncertainties of the predictions in Fig. 2 are 
from the NRQCD ME and PDF. The uncertainties from 
fi and rric are negligible. 

IV. SUMMARY 

In summary, we have presented the NRQCD predic- 
tions for the Pt and y distributions of the differential 
cross section and the polarization parameter a for the 
prompt J/ip produced in pp collisions at ^/s = 200 GeV. 
The Pt and y distributions are given in the kinematic 
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ranges 1.2 < |y| < 2.2 and 2 GeV < pr < 20GcV, re- 
spectively. While the NRQCD prediction for the prompt 
J/4' polarization in the midrapidity region {\y\ < 0.35) is 
slightly longitudinal, that in the forward-rapidity region 
(1.2 < \y\ < 2.2) is slightly transverse with less theo- 
retical uncertainties. The change in the lower-pr cut to 
2 GeV also reduced the uncertainties in our prediction 
by avoiding a possible break down of fixed-order calcula- 
tions at low pt- The LO CSM prediction severely under- 
estimates the differential cross section and shows more 
transverse J/tp than that of NRQCD. Our polarization 
prediction is less transverse as that of the s-channel-cut 
method [36] {a = 0.36), but is compatible with the NLO 
CSM prediction [32] , which has larger uncertainties than 
our results. We anticipate our NRQCD predictions pre- 
sented here can be tested against a forthcoming updated 



analysis by the PHENIX Collaboration. 
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